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Abstract 

An RPC carpet covering ~ 10 4 m 2 (ARGO-YBJ experiment) will be installed in the YangBaJing Laboratory 
(Tibet, China) at an altitude of 4300 m a.s.l. . A test-module of ~ 50 m 2 has been put in operation in this 
laboratory and about 10 6 air shower events have been collected. The carpet capability of reconstructing the 
shower features is presented. 
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1 Introduction 

As an Italian-Chinese collaboration project, the experiment ARGO-YBJ (Astrophysical Radiation with 
Ground-based Observatory at YangBaJing) is under way over the next few years at Yangbajing High Altitude 
Cosmic Ray Laboratory (4300 m a.s.l., 606 g/cm 2 ), 90 km North to Lhasa (Tibet, P.R. China). The aim of the 
experiment is the study of cosmic rays, mainly cosmic 7-radiation, at an energy threshold of ~ 100 GeV, by 
means of the detection of small size air showers at high altitude. 

The apparatus consists of a full coverage detector of dimension ~ 71 x 74 m 2 realized with a single layer 
of Resistive Plate Counters (RPCs). The area surrounding the central detector core, up to ~ 100 x 100 m 2 , 
consists in a guard ring partially (~ 50 %) instrumented with RPCs. These outer detector improves the 
apparatus performance, enlarging the fiducial area, for the detection of showers with the core outside the full 
coverage carpet. A lead converter 0.5 cm thick will cover uniformly the RPC plane in order to increase the 
number of charged particles by conversion of shower photons and to reduce the time spread of the shower 
front. The site coordinates (longitude 90° 31' 50" E, latitude 30° 06' 38" N) permits the monitoring of the 
Northern hemisphere in the declination band —10° < 5 < 70°. 
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ARGO-YBJ will image with high sensitivity atmospheric showers induced by photons with E~ t > 100 

GeV, allowing to bridge the GeV and TeV energy regions and to face a wide range of fundamental issues in 

Cosmic Ray and Astroparticle Physics (Abbrescia et al. (1996)): 
£5 , 

1) Gamma-Ray Astronomy , at a ~ 100 GeV threshold energy. Several galactic and extragalac- 
tic point candidate sources can be continuously monitored, with a sensitivity to unidentified 
sources better than 10% of the Crab flux. 

2) Diffuse Gamma-Rays from the Galactic plane, molecular clouds and SuperNova Remnants at 

energies > 100 GeV. 

3) Gamma Ray Burst physics , by allowing the extension of the satellite measurements over the 
full GeV/TeV energy range. 

4) p/p ratio at energies 300 GeV -j- TeV not accessible to satellites, with a sensitivity adequate 
to distinguish between models of galactic or extragalactic p origin. 

5) Sun and Heliosphere physics, including cosmic ray modulations at 10 GeV threshold energy, 
the continuous monitoring of the large scale structure of the interplanetary magnetic field and 
high energy gamma and neutron flares from the Sun. 

Additional items come from using ARGO-YBJ as a traditional EAS array covering the full energy range 
from 10 11 to 10 16 eV. Since the detector provides a high granularity space-time picture of the shower front, 
detailed study of shower properties as, for instance, multicore events, time and lateral distributions of EAS 
particles, multifractal structure of particle densities near the core, can be performed with unprecedented reso- 
lution. Detector assembling will start late in 2000 and data taking with the first ~ 750 m 2 of RPCs in 2001. 



In order to investigate both the RPCs performance at 4300 m a.s.l. and the capability of the detector 
to sample the shower front of atmospheric cascades, during 1998 for the first time a full coverage carpet 
of ~ 50 m 2 has been put in operation in the Yangbajing Laboratory. The present paper is a report of this 
test-experiment. 

2 The test experiment at Yangbajing 

The basic elements of ARGO-YBJ detector are RPCs of dimensions 280 x 125 cm 2 . The detector is 
organized in modules of 12 chambers whose dimensions are 5.7 x 7.9 m 2 . This group of RPCs represent a 
logical subdivision (Cluster) of the apparatus: the detector consists of 1 17 Clusters in the central part and 28 
Cluster in the guard ring for a total of 1740 RPCs. The proposed lay-out allows to achieve an active area of 
~ 92% the total. The trigger and the DAQ systems are built following a two level architecture. The signals 
from the Cluster are managed by Local Stations. The information from each Local Station is collected and 
elaborated in the Central Station. According to this logic a module of 12 RPCs (i.e. the Cluster) represents the 
basic detection unit. 

A Cluster prototype, very similar to that proposed for ARGO-YBJ experiment, has been installed in the 
Yangbajing Laboratory. It consists of 15 chambers distributed in 5 columns with an active area of ~ 90%. The 
total area is about 6.1 x 8.7 m 2 . 

The detector consists of single-gap RPCs made of bakelite (with volume resistivity p > 5 • 10 11 Q • cm) 
with a 280 x 112 cm 2 area. The RPCs read-out is performed by means of Al strips 3.3 cm wide and 56 cm 
long at the edge of which the front-end electronics is connected. The FAST-OR of 16 strips defines a logical 
unit called pad: ten pads (56 x 56 cm 2 ) cover each chamber. The FAST-OR signal from each pad is sent, via 
coaxial cable, to dedicated modules that generate the trigger and the STOP signals to the TDCs. Each channel 
of the TDCs measures, with 1 ns clock, the arrival times (up to 16 hits per channel) of the particles hitting 
a pad, for a total number of 150 timing channels. The RPCs have been operated in streamer mode with a 
gas mixture of argon (15%), isobutane (10%) and tetrafluoroethane C2H2F4 (75%), at a voltage of 7400 V, 
about 500 V above the plateau knee. The efficiency of the detector, as measured by a small telescope selecting 
a cosmic ray beam, is > 95%, and the intrinsic time resolution at ~ 1 ns. The description of the results 
concerning the RPCs performance at YangBaJing are given in Bacci et al. (1999). 

3 Data Analysis 

Different triggers based on pad multiplicity have been used to collect ~ 10 6 shower events with 0.5 cm of 
lead on the whole carpet in April-May 1998. The integral rate as a function of the pad multiplicity is shown in 
Fig. 1 for showers before and after the lead was installed. A comparison at fixed rate indicates an increasing 
of pad multiplicity due to the effect of the lead of ~ 15 -j- 20%, as expected according to our simulations. 

3.1 Detector Calibration The relative time offset among different pads are measured as follows: 

1) We construct the time distributions of each TDC channel adding all the delays in the individuals 

events and compute the relative time mean values tj. 

- y N _ u 

2) We compute the mean value of the U distribution: < t >= ^ - ■ 

3) We define the time offset as Atj = tl— < t >. These values are used to correct the times 

provided by each TDC channel. 

To check the consistency of the procedure, we fit on an event-by-event basis the corrected time values to a 
plane, construct for each TDC channel the distribution of time residuals Sti = t p i ane — U and calculate the 
mean values < Sti >■ These are distributed with a spread of ~ 0.6 ns. 

3.2 Event Reconstruction In this test-experiment we are not able to determine the core position, 
therefore we use a plane as the fitting function. Moreover, the estimated arrival direction is relatively free 
from the curvature effect because we sample only a small portion of the shower front. In this approximation 
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Figure 1: The integral rate as a function of the pad 
multiplicity. 



Figure 2: Time profile observed in a typical event. 
Straight lines are fit to experimental hits. 



the expected particle arrival time is a linear function of the position. The time profile observed in a typical 
event is shown in Fig. 2. Here x, y are orthogonal coordinates which identify the pad position. Straight lines 
are one-dimensional fits to experimental hits along two different x values. 

We performe an optimized reconstruction procedure of the shower direction as follows: 

1) Unweighted plane fit to hits for each event with pad multiplicity > 25. 

2) Rejection of out-lying points by means of a 2.5 a cut and iteration of the fit until this condition 

is not verified. 

After these iterations a fraction < 10% of the time signals that deviate most from the fitted plane are excluded 
from further analysis. The distribution of time residuals 5t = t p i ane — ti (Fig. 3) exhibits a long tail due to time 
fluctuations and to the curved profile of the shower front, more pronounced for low multiplicity. The width of 
these distributions is related to the time thickness of the shower front. Since the position of the shower core 
is not reconstructed, the experimental result concerns a time thickness averaged on different radial distances. 
Increasing pad multiplicity select showers with core near the detector, as confirmed by MC simulations. Taking 
into account the total detector jitter of 1.3 ns (RPC intrinsic jitter, strip length, electronics time resolution) the 
time jitter of the earliest particles in high multiplicity events (116 -r- 120 hits) is estimated ~ 1.1 ns. 

3.3 Angular Resolution The angular resolution of the carpet has been estimated by dividing the detec- 
tor into two independent sub-arrays and comparing the two reconstructed shower directions. Events with N 
total hits have been selected according to the constraint N oc id — N even ~ N/2. The even-odd angle difference 
A9 eo is shown in Fig. 4 for events in different multiplicity ranges. We note that these distributions narrow, as 
expected, with the increase of the shower size. To see the dependence of the angular resolution on the lead 
sheet, we show in Fig. 5 the median M/±o eo of the distribution of A9 eo as a function of pad multiplicity for 
showers reconstructed before and after the lead was added. The improvement of the angular resolution is a 
factor ~ 1 .4 for N = 50 and decreases with increasing multiplicity. 

Assuming that the angular resolution for the entire array is Gaussian (Alexandreas et al. (1992)), a standard 
deviation ag ~ 2.1° is found for events with a pad multiplicity > 100. 
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Figure 3: Distribution of time residuals for events 
with different pad multiplicity (all channel added). 



Figure 4: Even-odd angle difference distribution for 
events with different pad multiplicity. 



4 Conclusions 

A Resistive Plate Counters carpet of 



50 m 2 has been put in operation at the Yangbajing Laboratory to 
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study the high altitude performance of RPCs and the detector 
capability of imaging with high granularity a small portion of 
the EAS disc, in view of an enlarged use in Tibet (ARGO-YBJ 
experiment). 

In this paper we have presented the results of this test exper- 
iment concerning the carpet capability of reconstructing the 
shower features. In particular, we have focused on the angular 
resolution in determining the arrival direction of air showers, 
the most important parameter for 7-ray astronomy studies. 
The effect of a 0.5 cm lead sheet on the whole carpet has been 
investigated. An increase 15 -r- 20% of the hit multiplicity is 
found. The improvement of the angular resolution depends on 
the shower density. 

The test confirms that RPCs can be operated efficiently to sam- 
ple air showers at high altitude with excellent space and time 

resolutions. The results are consistent with data assumed in the figure 5: Median of At) 
computation of the performance of the ARGO-YBJ detector, function of pad multiplicity. 
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